Effects of formation properties in one-proton radioactivity 
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It is shown that the proton formation probability, extracted from experimental data corresponding 
to one-proton radioactivity, is divided into two regions when plotted as a function of an universal 
parameter. This parameter is derived from a microscopic description of the decay process. In 
this way we explain the systematics of proton emission half-lives. At the same time the formation 
probability is shown to be a useful quantity to determine the deformation property of the mother 
nucleus. 
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Understanding how nuclear many-body systems can 
self-organize in simple and regular patterns is a long- 
standing challenge in modern physics. The first case 
where this was realized was in a radioactivity, where the 
regularity manifests itself as striking linear correlations 
between the logarithm of the decay half-life and the ki- 
netic energy of the outgoing particle. This is known as 
the Geiger-Nuttall law, which was proposed 100 years ago 
[l[ . This law has been so successful that even today it is 
applied in studies of radioactive decays (see, e.g., Refs. 
@, yj)- The reason for this success is that the cv-particle 
formation probability, which is neglected in the Geiger- 
Nuttall law, usually varies from nucleus to nucleus much 
less than the penetrability. In the logarithm scale of the 
Geiger-Nuttall law the differences in the formation prob- 
abilities are usually small fluctuations along the straight 
lines corresponding to the isotopic chains. The impor- 
tance of a proper treatment of a decay was attested by 
a recent calculation which shows that the different lines 
can be merged in a single one Q . The resulting universal 
decay law (UDL), which uses three free parameters only, 
explains well all known ground-state to ground-state ra- 
dioactive decays. A similar three-term formula is also 
proposed in Refs. pl|(|. 

This good agreement is a consequence of the smooth 
transition in the nuclear structure that is often found 
when going from a nucleus to its neighboring nuclei. This 
is also the reason why the BCS approximation works so 
well in many nuclear regions. Notable discrepancies are 
only seen in a few cases around the N — 126 shell closure 
where clustering induced by the pairing mode is inhibited 
• Nuclei can even undergo more complex heavy-cluster 
decays, but these also follow similar linear relations 0- 
ill- 

Besides these cluster decay cases, proton radioactivity 
also provides a unique opportunity to study the structure 
of nuclei close to the proton drip- line [1(|. In the past 



decades tens of proton decay events have been observed 
in odd-Z elements between Z = 53 and Z — 83, leading 
to an almost complete identification of the edge of nu- 
clear stability in this region [ill, IHI ■ On the theoretical 
side, one may extract from proton emission detailed in- 
formation on the nuclear potential beyond the drip-line 
[I3I as well as the spectroscopic properties of the nu- 
clei involved fl5l - fT7| . 

The proton-emission process can be looked as a quan- 
tum tunneling through the Coulomb and centrifugal bar- 
riers of a quasistationary state jl8|. Similar to a and 
heavy cluster decays, the proton decay process can be 
divided into an "internal region" , where the compound 
state is restricted, and the complementary "external re- 
gion" . This division is such that in the external region 
only the Coulomb and centrifugal forces are important 
and the decaying system behaves like a two-particle sys- 
tem. The corresponding expression for the half life can 
be written as 
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where v and I are the outgoing velocity and the angular 
momentum carried by the outgoing proton, respectively. 
[L is the reduced mass corresponding to the final system. 
R is the radius dividing the internal and external regions. 
The half life does not depend upon R [l~8[ . At this point 
the wave function of the proton is matched with the cor- 
responding outgoing wave function in the external re- 
gion. In the following, the distance R will be taken as the 

1 /3 

touching point, i.e., R — Ro{AJ + 1), with i? =l-2 fm. 
The other quantities are standard, i.e., H^~(x,p) is the 
Coulomb- Hankel function with arguments \ — 2Zde 2 jhv 
and p = fivR/h. 

The amplitude of the wave function in the internal re- 
gion is the formation amplitude, 
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TABLE I. Coefficient sets of Eq. @ that determined by fit- 
ting to experiments of all proton decays (I) and to decays 
of nuclei with N < 75 and N > 75 separately (II), and the 
corresponding rms deviations. 





Emitter 


a 


b 


c 


d 


a 


I 


all 


0.386 


-0.502 


2.386 


-17.8 


0.440 


II 


N < 75 


0.443 


-0.364 


2.66 


-23.6 


0.344 


N > 75 


0.403 


-0.110 


2.766 


-27.8 


0.235 



where d, p and m label the daughter, proton and mother 
nuclei, respectively. VP are the intrinsic wave functions 
and £ the corresponding intrinsic coordinates. One sees 
from Eq. J2]) that J-i(R) would indeed be the wave func- 
tion of the outgoing particle ip p (R) if the mother nucleus 
would behave at the point R as 
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By applying a similar technique as in Ref. [4] and con- 
sidering the influence of the centrifugal barrier, one finds 
that the logarithm of the decay half-life can be approxi- 
mated by 



!)/(/ + c, (4) 
= A^Z p Z d Q- 1/ \ 



logT 1/2 = a X ' + bp' + dl(l-{ 

where a, 6, c and d are constants, x' 

p' = ^AZ p Z d {A d /3 + A 1 / 3 ), A = A d A p /(A d + A p ), and 
in this case of proton decay it is Z p = A p =1 (for details 
see Ref. Notice that c depends on the formation 

amplitude J-i(R). This formula should be compared with 
that of Ref. 17| where a linear correlation between the 
reduced half-life and the Coulomb parameter was found. 

The coefficients a to d can be determined by fitting 
available experimental data. For this purpose we took a 
total number of 44 decay events from the compilation of 
Ref. [H and the recent results of Ref. [2(J. The fitted 
values are listed in Table |TJ It is seen that the log values 
of the experimental half-lives can be reproduced within 
an error of a = 0.44. Correspondingly, the decay half- 
lives are reproduced within an average factor of three. 
This should be compared with the case of a decay where 
the UDL reproduces the available experimental half-lives 
within a factor of about 2.2. 

To probe Eq. (UJ), we plotted in the upper panel of Fig. 
[T]the quantity log T 1 J^' t ' — [bp' +dl(l+ l)/p'+ c] as a func- 
tion of x' ■ I n the lower part of this Figure we plotted the 
discrepancy between experimental and calculated values, 



i.e., the ratio 1Z = T, 



Expt. 
1/2 



as a function of the emit- 



ter charge numbers Z . The calculations were performed 
by using the parameter set I in Table HI It is seen that 
most of the data can be reproduced by the calculation 
within a factor of 4, i.e., with 0.25 < 1Z < 4. Larger dis- 
crepancies are seen for a emitters between 63 < Z < 67 
and the isomeric /in/2 hole state in the Z — 81 nucleus 
177 T1, where the experimental decay half life is underes- 
timated by the calculation by a factor of about 8. 
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FIG. 1. (Color online) Upper: UDL proton decay half-lives. 
Dots are experimental values while the straight line is ax' 
with a provided by the parameter set I of Table I. e is given 



Expt. /rpCal. 
1/2 / 1 l/2 



by e = bp' + dl(l + l)/p' + c. Lower: 11 = T : 
a function of the charge number Z. Within the dashed lines 
is the region where the data is reproduced by the calculation 
within a factor of 4. 



A better understanding of these discrepancies would 
require a systematic study of the formation probabilities. 
We performed this task by extracting the formation am- 
plitudes F(R) from the experimental half lives by using 
the expression Q, 



log\RF(R)\- 2 = logT^f- 



- log 



— \Ht(x,P)\ 



(5) 

which we plotted in Fig. [2] as a function of p' . 

One sees in this Figure that two clearly defined regions 
emerged, as indicated by the dashed lines which we ob- 
tained by using the parameter set II of Table Q] 

The origin of these two regions can be gathered from 
the values of the Q values, deformation parameters and 
formation probabilities presented in Table |nl where the 
experimental and corresponding theoretical half-lives are 
also given. 

One sees in this Table that the region to the left in 
Fig. [21 i.e., for lighter isotopes, corresponds to the de- 
cays of well deformed nuclei. The formation probabilities 
decreases for these nuclei as p' increases. Then, suddenly, 
a strong transition occurs at /o'=20.5, corresponding to 
the nucleus 1 ||Tm. Here the formation probability ac- 
quires its maximum value, and then decreases again as 
p' increases. One notices in Table |TT] that this tendency 
is followed by the deformations, which gradually dimin- 
ishes after the transition point. One thus sees that the 



3 



TABLE II. Experimental and calculated half-lives for proton decays to ground states. Experimental Q values and half-lives 
are taken from Refs. [13, H3|. Stars in the emitters indicate isomeric states. The deformation parameter ft are from Ref. [2lj . 
The calculations were performed by using the parameter sets I (Ty£'(I)) and II (Tw^'(II)) of Table U In the last two columns 
are the formation amplitudes \Fi(R)\ and formation probabilities \RFi(R)\ 2 (Eq. (HI). Units in the calculated half lives are as 
in the corresponding experimental values. 



Emitter 


I 


Jm 


Q(keV) 


ft 


TiExpt. 

1 1/2 


T l/2 (I) 


T-iCal. / TT\ 

T l/2 ( n ) 


|F((i?)|(fm °i ) 


\RFi(R)\ (fm ) 


'sll 


2 


5/2+ 


827( 5) 


0.160 


93.5(5) fis 


243 


187 


0.020 


0.018 


l 5lCs 


(2) 


(0+3+) 


823( 7) 


0.208 


0.5(1) ms 


1.18 


1.39 


0.021 


0.021 




2 


3/2+ 


976( 3) 


0.207 


16.7(7) /us 


15.5 


9.72 


0.010 


0.005 


l 57La 


2 


3/2+ 


814(11) 


0.290 


23.8(20) ms 


6.52 


12.3 


0.008 


0.003 


l 57La 


4 


9/2+ 


951( 6) 


0.290 


10(5) ms 


7.55 


12.6 


0.016 


0.013 


mPt 


2 


3/2+ 


900(10) 


0.318 


10(H) ms 


1.72 


3.28 


0.006 


0.002 


mEu 


2 


1+ 


1039(15) 


0.331 


0.90(1^) ms 


0.474 


1.12 


0.010 


0.005 


kjEu 


2 


3/2+ 


959( 9) 


0.331 


21.4(1^) ms 


4.18 


13.8 


0.007 


0.002 


MfTb 


3 


7/2" 


1200( 7) 


0.325 


0.94(1°;") ms 


0.181 


0.428 


0.006 


0.002 


x #Ho 


3 


(6", 0", 8+) 


1106(10) 


0.297 


6(3) ms 


4.91 


23.3 


0.016 


0.014 


wHo 


3 


7/2~ 


1190( 8) 


0.286 


4.1(1) ms 


0.650 


2.30 


0.006 


0.002 


^Ho* 





1/2+ 


1255( 8) 


0.286 


6.6(1°,;?) ^s 


6.04 


12.0 


0.010 


0.006 


raTm 


5 


(10+, 5") 


1725(16) 


0.258 


2.7(1J;?) MS 


22.0 


6.37 


0.045 


0.111 


mTm 


5 


11/2" 


1753( 7) 


0.249 


~3.46(32) /is 


14.8 


4.28 


0.031 


0.056 


^Tm 


5 


(5") 


1210( 4) 


-0.199 


117.6(64) ms 


151 


67.8 


0.030 


0.051 


"fTm* 


5 


(10+) 


1140( 4) 


-0.199 


203(6) ms 


792 


383 


0.058 


0.190 




5 


11/2" 


1073( 5) 


-0.190 


3.78(1.27) s 


4.36 


2.32 


0.035 


0.069 


1 A7m * 


2 


3/2+ 


1133( 3) 


-0.190 


0.360(36) ms 


1.51 


0.269 


0.031 


0.056 


"iLu 


5 


> 5 _ 


1283( 3) 


-0.164 


64.0(56) ms 


83.1 


51.3 


0.030 


0.050 


x ??Lu* 


2 


(1-.2-) 


1306( 5) 


-0.164 


43(11) ^ 


94.1 


20.9 


0.020 


0.023 


n'Lu 


5 


11/2" 


1253( 3) 


-0.156 


127.1(18) ms 


155 


100 


0.030 


0.051 


^Lu* 


2 


3/2+ 


1332(10) 


-0.156 


16(1) fis 


54.0 


11.9 


0.024 


0.033 




5 


11/2" 


1468(15) 


0.008 


2.9(1^) ms 


5.56 


4.37 


0.031 


0.056 


^fTa 


2 


(2~) 


1032( 5) 


-0.050 


149(8) ms 


286 


135 


0.029 


0.050 


^fTa* 


5 


(9+) 


1127( 7) 


-0.050 


8.52(2.12) s 


9.19 


10.3 


0.036 


0.076 


"jTa 





1/2+ 


947( 7) 


0.045 


0.300(105) s 


0.941 


0.401 


0.037 


0.080 


'IfRe* 


5 


11/2" 


1831(20) 


0.053 


20.2(37) (j>s 


47.7 


42.6 


0.026 


0.041 


're Re 


2 


(2-) 


1287( 6) 


0.080 


0.687(11) ms 


1.22 


0.642 


0.021 


0.027 


re 1 Re 





1/2+ 


1214( 6) 


0.080 


0.440(2) ms 


1.45 


0.660 


0.027 


0.045 


'Is Re* 


5 


11/2" 


1338( 6) 


0.080 


224(31) ms 


169 


226 


0.024 


0.034 


164t 
77 lr 


5 


(9+) 


1844( 9) 


0.089 


0.113(l|g) ms 


0.0829 


0.108 


0.015 


0.014 


165t * 
77 lr 


5 


11/2" 


1733( 7) 


0.099 


0.34(7) ms 


0.385 


0.550 


0.021 


0.026 


166t 
77 ir 


2 




1168( 7) 


0.107 


0.152(71) s 


0.0678 


0.0631 


0.014 


0.011 


166t * 
77 lr 


5 


(9+) 


1340( 8) 


0.107 


0.84(28) s 


0.419 


0.835 


0.021 


0.026 


^Ir 

77 11 


o 


1 /2+ 

1/ z 




0.116 


1 1 0(1 5") ms 

-L '.J I Illo 


112 


0.919 


0.020 


0.023 


^Ir* 


5 


11/2" 


1261( 7) 


0.116 


7.5(24) s 


2.45 


5.39 


0.019 


0.021 


^Au 


2 


(2-) 


1488(12) 


0.080 


321(1{£) /xs 


153 


1.53 


0.010 


0.007 




5 


(9+) 


1770( 6) 


0.080 


1.046(11^) ms 


0.482 


0.991 


0.014 


0.011 


V^Au 





1/2+ 


1464(10) 


-0.105 


24.5(1 4 ;I) ^ 


54.8 


45.1 


0.020 


0.025 


^Au* 


5 


11/2" 


1719( 4) 


-0.105 


2.22(19) ms 


1.01 


2.18 


0.014 


0.012 


81 11 





(3", 4", 5") 


1282(18) 


-0.053 


5.2(1?;°) ms 


8.55 


12.7 


0.023 


0.032 


177 T i 
81 11 





1/2+ 


1180(20) 


-0.053 


67(37) ms 


119 


203 


0.027 


0.045 


177 T] * 
81 11 


5 


11/2" 


1984( 8) 


-0.053 


396(1^) fj>s 


50.1 


137 


0.006 


0.003 


185 R - 





1/2+ 


1624(16) 


-0.052 


58(9) (is 


19.0 


33.1 


0.008 


0.004 
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FIG. 2. (Color online) Proton-decay formation amplitudes 
log 10 \RF(R)\ 2 extracted from experimental data as a func- 
tion of p . Squares correspond to nuclei with N < 75(Z < 67) 
while circles are for N > 75(Z > 67). 
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FIG. 3. The parameter p versus the fragmentation potential 
(Eq. ([6])) for nuclei with Z < 67 (open circles) and Z > 67 
(dark circles). 



reason of this tendency of the formation amplitude is 
related to the deformation. Indeed, in the left region 
of Fig. [21 the decays of the deformed nuclei proceed 
through small spherical components of the corresponding 
deformed orbitals and, therefore, the formation probabil- 
ities are small. 

The right region of Fig. [2] involves the decays of spher- 
ical orbits as well as major spherical components of de- 
formed orbitals (for example, ft.xi/2 component of the or- 
bital 11/2 [505]). 

Another striking feature is that the formation proba- 
bility of the odd-odd nucleus 144 Tm is abnormally large. 
Considering that the experimental uncertainties regard- 
ing the half-life (from where the formation probabil- 
ity is extracted) is large one may doubt the validity of 
this anomaly. Moreover, it is not surprising, given this 
anomaly, that for 144 Tm the half-life evaluated with the 
set of parameters I differs from the corresponding exper- 
imental value by a factor R which is beyond the scale 
used in the lower panel of Fig. [1] 
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FIG. 4. (Color online) The formation amplitudes |fi(i?)| ex- 
tracted from experimental data for proton decays of nuclei 
N > 75(Z > 67) as a function of u calculated from BCS cal- 
culations using for the Woods-Saxon mean field the universal 
parameters [23[ (upper) and the Cherpunov parameters [24[ 
(lower). 



It is interesting to point out that a similar "clustering" 
of the proton wave functions into two distinct regions, 
as the one seen in Fig. [21 was evidenced in Ref. (22j. 
But there the determining variable was not p' but the 
so-called fragmentation potential. This is given by the 
difference between the Coulomb barrier and the Q-value 
as 



Vfrag = 



Z d e 2 
R 



(0) 



This feature is not surprising, since there is a linear cor- 
relation between these parameters, as can be seen in Fig. 
[3] where the parameter p' is plotted versus the above de- 
fined fragmentation potential. One sees here a clear lin- 
ear correlation between these parameters. Moreover, the 
plot is divided into two regions with Z < 67 and Z > 67. 
In fact that correlation is stronger than what the figure 
suggests since the upper part of the plot (dark circles) is 
divided into several regions where few close values of the 
fragmentation potential (with different Q values) corre- 
spond to the same p' . It is worthwhile to point out that 
the linear correlation between Vf rag and p' , as seen in Fig. 
[3l is a result of the matching of the proton internal wave 
function with the corresponding external Coulomb wave. 
This is because both parameters are connected by this 
procedure. 

Going back to Fig. [5] we notice that in the BCS 
approach the formation amplitude at a given radius R 
would be proportional to the product of the occupancy u 
times the single-proton wave function tp p (R). Therefore 
the tendencies seen in Fig. [5] may be due to the BCS am- 
plitudes or the radial wave functions. To recognize the 
influence of the value of the BCS amplitudes we plotted in 
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FIG. 5. (Color online) Single-proton wave functions in Lu 
for different channels calculated with the universal Woods- 
Saxon parameter. The inserted plot shows the wave func- 
tions of the 0/in/2 orbital with different energies derived by 
changing the depth of the potential. 



Fig. |4]the formation probabilities extracted from 

experiment for the case of proton decays corresponding 
to nuclei with N > 75(Z > 67) as a function of u. The 
u values were calculated by using a Woods-Saxon poten- 
tial according to the universal (upper) and Cherpunov 
(lower) parameters. One sees that the tendency is the 
same for both sets of parameters, namely that the val- 
ues of of the formation probabilities increases with it, 
strengthening our interpretation of the behavior of the 
formation probabilities presented above. 

Still one has to consider eventual changes in the radial 
wave functions at the nuclear radius (where the formation 
probabilities are determined) . For this we plotted in Fig. 
[S]the single-particle wave functions <fr(r) = rip p (r) for dif- 
ferent orbitals corresponding to the nucleus 151 Lu. It is 
seen that the values of <fi(r) around the nuclear surface 
are quite similar to each other. From this we conclude 
that the fluctuations in the experimental formation am- 
plitudes found above for nuclei with N > 75 are mainly 
due to fluctuations of the it values. It is also to be noticed 
that some cases which departs from the UDL correspond 
to the decays of hole states. This occurs in the isomeric 
state of 177 T1 (as already pointed out above) and also in 
the ground state of 185 Bi [25j . 



Another feature that has to be considered is the influ- 
ence (if any) of the Q value upon the formation proba- 
bilities. The Q value determines the penetrability and, 
therefore, the radioactive decay process. The question 
is whether even the spectroscopic quantities are affected 
by the Q value. In Ref. [22| it was evidenced the fact 
that the logarithm of the wave function is proportional 
to the fragmentation potential (Eq. for all decay 

processes, in particular for proton emission. This is a di- 
rect consequence of the one-body Schrodinger equation. 
Anyway, it seems that the dependence on the Q-value 
alone is small in proton decay processes. To analyze this 
we show in the inserted plot of Fig. [5] the wave functions 
of the Qhxx/2 orbital under different energies. These were 
obtained by changing the depth of the potential. As per- 
haps expected, the formation amplitude at the nuclear 
surface is not sensitive to changes in the energy, i.e., in 
the Q values. Neither the amplitudes u are much affected 
by the changing of the potential depth, as also expected. 
One therefore concludes that the changes in the forma- 
tion amplitudes extracted from experimental data are a 
result of changes in the deformation of the decaying nu- 
clei. 

In conclusion, we have shown in this paper that the 
proton formation probability in the corresponding decay 
process depends upon the deformation of the decaying 
nucleus. In a well deformed nucleus the decay proceeds 
through one of the spherical components of the deformed 
orbit, which is usually small in this case of large deforma- 
tions. Therefore the formation probability is small. On 
the contrary, in spherical or weakly deformed nuclei the 
decay proceeds through the only component that is avail- 
able and, as a result, the formation probability is large. 
We have strengthen this interpretation by showing that 
in the case of weakly deformed nuclei the formation prob- 
ability is proportional to the BCS occupation number u 
corresponding to the decaying orbit. 
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